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Abstract 

Local and global existence of localized solutions of a discrete non- 
linear Schrodinger (DNLS) equation, with arbitrary on-site nonlinear- 
ity, is proved. In particular, it is shown that an initially localized 
excitation persists localized during infinite time. Moreover, if initial 
localization is stronger than \n\~'^ with any power d, it maintains it- 
self as such during infinite time. The results are generalized to various 
types of inter-side and saturable nonlinearities, to lattices with long 
range interactions, as well as DNLS with dissipation. 

Key words: Discrete nonlinear Schrodinger equation; localized solutions, local and 
global existence. 



1 Introduction 



Localization and transfer of energy in nonlinear lattices are the key problems of the 
physics of discrete systems, which have attracted a great deal of attention during 
the last years (see e.g. jj). In this context, the discrete nonlinear Schrodinger 
(DNLS) equation 

iqnit) + A(Zn(t) - Xkn(t)p(^-')gn(t) = , (1) 
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where p is an integer which throughout this paper will be considered bigger than 
or equal to 2, p > 2, x = il, and Ag„ = Qn+i + Qn-i — '^Qn, is an object of special 
interest P] due to its numerous physical applications, ranging from dynamics of 
electrons in one-dimensional crystals in the tight-binding approximation, to arrays 
of optical waveguides, with Kerr (p = 2) and non-Kerr (p > 3) nonlinearities, and 
arrays of Bose-Einstein condensates. 

A possibility of energy localization in the form of so-called intrinsic localized 
modes, i.e. periodic in time and strongly localized in space solutions of Eq. 
was predicted in Ref. |3j. Later on, in Ref. this fact received its mathematical 
proof, and recently a classification of diversity of intrinsic localized modes has 
been provided 5 . By its definition, an intrinsic localized mode can be viewed as 
a particular solution of the Cauchy problem for the DNLS equation, which has 
sufficiently large amplitude and strong spatial localization, thus belonging to a 
subclass of all possible solutions. Meanwhile, the problem of localization properties 
of an arbitrary initially localized pulse, remains open. This justifies one of the main 
objectives of the present paper: to provide a rigorous proof of persistence in time of 
localization of solutions of the Cauchy problem for the DNLS equation subject 
to an initial condition 



where a„, with n G Z, are given complex constants. Initial localization of the pulse 
will be characterized by the constraint 



where d is a positive integer, d > 1, which characterizes the decay of \am\ when \m\ 
tends to infinity. If ^ is verified for all d < do and fails for d > do the localization 
is referred to as power. As it is clear, for modes localized exponentially (which 
is also the case of intrinsic localized modes mentioned above), © holds for all 
positive d. 

Another important issue related to the DNLS model is that it can be viewed 
as a discretization of the nonlinear Schrodinger equation with high nonlinearity. 
As it is known jo., such equation at a given sign of a possesses blowing up solutions. 
Collapse occurs due to small wave-lengths (large frequencies). Thus it is natural to 
expect that the discretization, which introduces the linear spectrum frequency cut- 
off, prevents collapse in with any integer p and both x = il, which indeed was 
observed in 0. In the present paper consideration of the localization properties of 
solutions of the DNLS equation will be reduced to the existence theorems and thus 
arresting collapse will be automatically guaranteed (we notice that this approach is 
analogous to one used in [S] for nonlinear Klein-Gordon and sine-Gordon lattices). 

The organization of the paper is as follows. In Sec.|2l we describe the statement 
of the problem, define an appropriate Banach space, and formulate the main result. 
Sec. 01 is devoted to the proof of the local existence of localized solutions. Global 



gn(0) = a, 





(3) 



m=—oo 
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existence is proved in Sec. \^ In Sec. [5] the results are generalized to the DNLS 
equation with dissipation, as well as models with long-range interactions and other 
types of nonlinearity. The results are summarized in Conclusion. 



2 Statement of the problem and main results 

By analogy with |S] our proof will be based on the Banach fixed point theorem jH] 
(see also below). For its implementation we rewrite the nonlinear model (pQ) in an 
integral form. To this end we first consider the linear problem associated with (^: 



ipn{t) + Apn{t) = , pn{0)=an , (4) 

which solution is trivially obtained in the form 

+00 

Pnit) = ^ amG{n-m,t), (5) 

m=— oo 

with G{n, t) being the Green function 



1 /■a 



Try.. 

Then, using © and © we rewrite as 

-oo „J 



qn{t)=Pn{t)+ ^ / G{n-m,t- s)F{qm{s))ds 

m=-oo 



(7) 



where F {qn{t)) is given by: 

Fiqn{t))=x\qn{t)\^^''-'\n{t) . (8) 

One of the main objects of our analysis will be a 2(i-momentum Qd{t) defined 
as follows 

+ 00 

QM= {l + n'Y\qn{t)\' . (9) 

n=— oo 

At (i = 0, © becomes a trivial integral of (Q), 



+00 



Q0= E kn(t)|^ ^ = 0, (10) 

n=— 00 

which often is referred to as a number of particles N and therefore below we use 
notation Qq = N. For an arbitrary d Qd{t) may grow, but it always remains 
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bounded. Indeed, one can use Cauchy-Schwarz's inequality and compute a deriva- 
tive of Qdit) from (P), in order to deduce the estimate as follows: 



d 
di 



+00 ^ d 



+ CXD 



<Co {^ + i^+nfy^\l + n^f^\qn\\qn+i\<CoQd{t). 

ra=— 00 

Here Co = [{Vb + if - (\/5 - lf]/2'^-'' and it was used 

V5-1 1 + {1 + nf V5 + 1 

Then the positive real function Qd{t) satisfies the hypothesis of Gronwall's 
lemma and hence 

Qd{t) < gd(0)e^o* . (11) 

The next step to be performed, for the sake of the use of the Banach fixed 
point theorem, is defining the appropriate Banach space. To this end, for some 
T > 0, we introduce a function q{t): 

q{t) = {..., qn-i{t),qn{t),qn+i{t),...) , (12) 

such that: qn{t) are continuously differentiable for any n, qn{0) = a„ and 

sup I 2 {l+n^y\qnit)A<^, 

and consider the linear space X'^{T), which consists of all such functions q{t). We 
also define a norm of an element q{t) € X'^(T) by 



hoo 



||g||^.= sup (1 + ^') kn(t)|' . (13) 

Then the space [X'^{T), \\-\\xd) becomes a Banach space. 

As it is evident, the norm can be used as a measure of localization of 

the solution: if a solution belongs to a Banach space X'^{T) with larger d, then its 
energy displays higher localization. 

We also observe that if 

Qd= sup ( V n^''\qn{t)A <oo , (14) 
0<t<T 
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for some integer d, then 

f +00 



sup 

0<t<T 



^ {1 + \qn{t)\H < 00, 



because Q^i < Qd for all di < d. This will allow us to simplify analysis by restrict- 
ing consideration to Qd whenever convergence of sums will be questioned. 
One can associate the map P to the integral form of Eq. ((71): 

+00 „t 

Pqn{t) = Pn{t) + Yl / G{n-m,t- s)F{qm{s))ds , (15) 



and respectively Pq{t) = . . , Pqn-i{t), Pqn{t), Pqn+i{t), . . .j . Then a solution 

of (P) is nothing but a fixed point of the map P: q{t) = Pq{t). This constitutes 
a basis for the use of the Banach fixed point theorem which states that: if one 
has a contraction mapping P from a closed subset F of a Banach space X'^ onto 
F, then there exists a unique q G F such that Pq = q. Thus, we have to define 
a closed subset of the Banach space as well as a functional that maps this subset 
into itself, the latter being a contraction. To this end we consider R > and p{t) 
as m (tra . and define the following subset of X'^{T): 

Yr{T) = {g G X'^(r), \\q - p\\^, < R, (/„(0) = ^,(0) = a„} . 

Clearly, Yr is a closed subset of X'^{T). Moreover, we observe that if g G Yji{T), 
then \\q\\xd < WpHxd + R < 00, and n^'^|(ir„(t)p < \\q\\xd < 00. 

Now we can formulate the main results of the paper, namely the existence 
theorems. 

Theorem 2.1. (Local existence) Assume that a„ satisfies the condition (0), then 
there exists T > and a unique function q{t), as in defined in [0, T[ which 

belongs to the Banach space X'^iT), such that qn{t) satisfies Q). 

Theorem 2.2. (Global existence) The local solution of ^ found in Theorem \2.1\ 
can be extended for any T > 0. Such global solution is unique. 



3 Proof of the local existence theorem 

In this section we prove Theorem 12.11 To shorten notations we suppress the 
subindex X'^ in the norm definition (i.e. || • \\xd will be designated as || • ||), 
taking into account that no other norm will be used. Also we denote constants 
(independent on the site number n and time t) by C, whose values will however 
vary from line to line. 

Let us consider some auxiliary results starting with 



5 



Lemma 3.1. For t G [0,r[ 

10(0,^)1^ = 1, and |G(n,t)|^<^ for n ^ 0, \/j e Z+ , (16) 

+ 00 

where Kj are constants depending on j and T , and the series n^'^|G(n, t)|^ 

n=—oo 

is uniformly convergent. 

Proof. For n = the estimate for the Green function trivially fohows from its 
definition Q. For n 7^ integrating by parts we obtain: 

G{n,t) = ^^ f" flj\a,t)da , 

''~2 

where fa\a,t) is the j*^ derivative of the function f{a,t), defined in ©, with 
respect to a. The derivatives can be represented as: 



fi^\a,t) = Y,M^)t''e 



Ait sin a+2ian 
h=l 



Here the functions (/>/i(ct), h = 1, . . . ,j are functions of sin(o") and cos((t), only, and 
thus (ph{cr) are periodic bounded functions of a. So we can deduce for t £ [0, T[: 



2 ^ , \ , / 1 ^ ^ 1 



2j 



where Ch = {(^hMlda. Now it is immediate to obtain (|T6l) and derive the 

71"^ 7-77/2 

1 ^ 

uniform convergence of the series: indeed, it is sufficient to take Kj = ChT^, 

h=l 

and j > d. □ 

It should to be mentioned here that the above result is nothing but a conse- 
quence of the Riemann-Lebesgue lemma (see e.g. ^Oj) stating that 

/(a,t)da = o(ni-^) , 

when n — > 00, i.e. that the integral decays faster that any power of n. Thus for all 

+00 

d > we always can find sufficiently large such that the series n'^^\G{n, i)|^ 

ri=— 00 

is convergent. 
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From Lemma l3. II we deduce a corollary important for the next consideration: 

oo 

^ n^''\G{n-m,t)\^ <Cd{l + mY , (17) 



where is a constant depending on the power d. 

For the next consideration it is also important that the localization principle 
holds also for the linear problem @, what is guaranteed by the following 



Lemma 3.2. The function: 

p{t) = (. . . ,Pn-l{t),Pn{t),Pn+l{t), . . .) , 



(18) 



where Pnit) is the solution for the linear problem belongs to the Banach space 
X'^iT), for any T > 0, provided that ^ holds. 

Proof. Using with j = d + 1, Cauchy-Schwarz inequality, and inequality (|17|) 
we obtain: 



,2d 



n 



n=—oo 



n=— oo 



amG{n-m,t) 



m=~oo 



-oo 



+ 00 

mi,m=— oo n=— oo 



n'^\G{n-mi,t)\-n'^\G{n-m,t)\ 



< 



oo / +00 

m=—oo \n=— oo 



J2 E n^^lGin-m^tf 



1/2- 



1 2 



+ 00 



< 



vm=— oo 



To conclude the proof we take into account that boundness of vl^ (see (jSJ) implies 
convergence of the last sum. □ 

Now we are in a position to present the proof of the main Theorem 12.11 

Proof of Theorem \2.1\ Let q{t) G Yr(T), and P is the map defined above. In order 
to prove that Pq{t) E Yii[T) we have to show that \\Pq — p\\ < R. Let 



+ 00 

(t) = Pqnit) - pn{t) = Y / Gin-m,t- s)F{q^{.s))ds 



m=—oo 



Thus we have to estimate (to shorten notations hereafter we use t = t — s and 
Ti = t — si and overbar hereafter stands for complex conjugation): 



hoo 



2d 



n 

n=—oo 

-oo „i 



/ G{n - m,T)\qm{s)\^'~^~^\m{s)ds 



n - mi,Ti) |gmi(si)p^^ 9mi('5i)c^si 
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We will verify later on that the series which appear above are uniformly con- 
vergent so that we can exchange order of summation and integration and estimate: 



< ds dsi V] r?'^ 

+00 

^Yy(lra{.s)\^^'^ \G{n - m,r)| ^ kmi (si)|^^"^ \G{n - mi,Ti) 



n=— 00 

+00 

X 

m=— 00 



mi=—oo 



r-t ft +0° 

Jo Jo ™_ „ „ 



2p-l 



10 Jo 



m=—oo 
-3 



To obtain the last inequality we used (|17|) . So we can conclude that 

\\wf < C\\qfN^P-^T^ < C {\\p\\ + Rf N^P-^T^ . 

Therefore we have only to choose T > sufficiently small, in such a way that: 

R 



T < 



C(||p|| +i?)^A^P-3/2 



(19) 



Then \\Pq — p\\ < R, that is the functional P maps Yji(T) into itself. 

Now we will prove that P is a contraction, that is, there exists a real number 
a, < a < 1, such that: 

\\Pq - Pq\\ < a\\q - q\\ Vg,gGYR(T). 

To this end we define a function 

Qmit) = Qmit) " k^Wl'^^"') qm{t) , 

and estimate 

hoo „i 



Y: n''\Pqn-Pcin\'= E 



n 



2d 



y ^ / G{n — m,T)Qrn{s)ds 
.—^Jo 



m=— 00 



+00 nt _ _ 

/ G{n- mi,Ti)Qrm{si)dsi 

,=-r^JO 



mi=— 00 



< ds dsi ^^'^ |G(n-m,r) 

•'0 ?i=-oo m=-oo 

+00 

X X] |0mi(si)| [^(n - mi,Ti)| . 
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Likewise in the previous part we have exchanged order of summation and integra- 
tion. 

Next we notice that 

2p-2 

\@nM < \Qm{t) - qm{t)\ ^ \qm{t)\^''~^~''\qm{tt , (20) 

fc=0 

and thus 

+00 



m=—oo 
+00 

< E (l + m^) \qni - qm\^ {\qm\ + \qm\)'^^~'^ 
m=— 00 



where N = l^mP- Thus we obtain the following estimate 

+00 +00 +00 

E n^^^ \eM\\G{n-m,r)\ \&m^{si)\\G{n - mun] 

n=—OD m=— 00 mi=— 00 

+00 +00 / +00 \ ^/"^ 

Y \@m{s)\\emAsi)\ I E ^''|G(n-m,r)n 



m=— 00 mi =—00 \n=— 00 

/ +00 \ 1/2 

|2 1 



X 



+00 +00 



J2 J2 a + mY^\i + mlf/^\em{s)\\emAsi)\ 

m=—oo mi=— 00 

<C\\q-qf{\\q\\ + \\q\\f , 
where has been used. From the last result we achieve: 

\\Pq-Pqf < sup { f ds f dsiC{\\q\\ + \\q\\)^\\q-q\ 
0<t<T Uo Jo 

< C{\\p\\+R)'\\q-qfT^ . 

Therefore given a positive number a, such that < a < 1 we can choose T 
small enough so that C {\\p\\ + R)'^ T < a, that is : 

T < 1 , (21) 

C{\\p\\+R)' 
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which imphes that P is a contraction in Y/j. Therefore choosing T = Tq that 
satisfies the conditions (fT9|) . (|2T|) we can conclude that there exists a unique fixed 
point q S Yji{T) of the map P, i. e. Pq = g for < t < Tq where q{t) is given by 
(fT2|l and g„(t) satisfies 0. □ 

One has to note that from the maj orations achieved during the proof one 
deduces that the series of functions are majorated by convergent numerical series. 
Then, from Weierstrass convergence criteria, one asserts that the series exhibit 
absolute convergence as well as uniform convergence, allowing one to change order 
of integration and summation. 

4 Global existence of solutions 

In order to prove global existence we have to extend the solution ((T)) for any T > 0. 

Proof. To this end we can use Zorn's lemma i.e. to extend the solution to 
the maximal interval of existence which we denote by [0, Tmaxi, and to prove that 
Tmax = +00. We proceed as follows: we suppose that Tmax < +oo and \\q\\ +oo 
as T — > Tmax and then we will get a contradiction by an a-priori estimate. This 
is a trivial task since from inequality we deduce: 

\\qf = sup Qd{t) < sup {Qrf(0)e^o*} = Qd(0)e^°^ < oo . (22) 

0<t<T 0<t<T 

Thus Tmax = +00 and the solution Q is a global solution of problem 0. 

It remains to prove the uniqueness of the global solution. Let rn(i) and 
be two solutions of the problem (0, @ and define Zn{t) = rn{t) — Sn{t). We will 
prove that Zn{t) = ,Vt > 0, what means that problem have a unique 
solution. To this end we consider the following function: 

Z{t)= {l + n^Y\z^{t)\\ (23) 

ra=— oo 

and compute its derivative in the form Z[t) = Di + D2 where 

{+00 
^ (1 + n^)'^ [zn+l{t)Znit) + Zn-l{t)Znit)] 
n=—oo 

{+00 
^ {l + n^yfn{t)Sn{t) (|r„(i)p(^-^) - |Sn(t)p(^-')) 
n=— 00 

With a similar argument used to estimate Qdit) we obtain Di < CoZ{t). Next, 
using that |9(f„s„)| < |r„| |r„ — Sn\, |?^nP < -^j and < A^, we estimate 

+00 

D2 < 2 {l + n^y\Q{r^Sn)\\rn-Sn\{\rn\^P-'' + --- + \snf^-') 

n=~oo 

< 2{2p - 2)N'^P-'^Z{t) . 
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Therefore Z{t) satisfies tlie following relation: 

Z{t) < CoZ{t) - A{p - l)N^P-^Z{t) = CZ{t) , 

and due to Gronwall's Lemma we have: 

< Z{t) < Z(0)e^* . 

Since Z{0) = (r„(0) = s„(0) = an) we thus have that Z(t) = ,\/t > and 
accordingly = , Vt > 0. □ 

5 Generalizations 

5.1 Dissipation 

The results obtained above are directly generalized to the DNLS equation with 
dissipation 

iilnit) + AUn{t) - x\Un{t)\'^^''~'^\n{t) + i-/Un{t) = , U„(0) = a„ . (24) 

where 7 is a positive constant characterizing dissipation. To this end we notice 
that by the means of the ansatz M„(t) = e~'^^qn{t) the problem H24() can be reduced 
to: 

iUt) + Ag„(t) - xe~^^''''^^'\qn{t)?^''-'\n{t) = , (?n(0) = a, . (25) 

The only change of the proofs presented above are in the definitions of F (see (jSJ) 
which now reads 

F {quit)) = e-^^P-'^^'x\qnit)\^^^-'Un{t) . (26) 

It is also important that in the case 7 > one has N(t) = e~^^N, where N is like 
in and thus an estimate < N. Thus, without presenting a proof, we 

formulate the following theorem 

Theorem 5.1. Assume that an, n G satisfy the condition then there exists 
T > and a unique function u{t) = {. . . ,Un-i{t),Un{t),Un+i{t), . . .), defined in 
[0,T[ which belongs to the Banach space X'^{T), such that Un{t) satisy 0). This 
solution ^ can be extended for any T > 0. Such global solution is unique. 

5.2 Long range interactions 

The results expressed in Theorems 12.11 and 12 . 21 are essentially based on the proper- 
ties of the Green function expressed by the formulas ((TB|) and (fTTI) . Let us consider 
now a generalized discrete nonlinear Schrodinger equation 

00 

iqnit)+ Yl K{n-m)q„,{t)-x\qnit)\^^'"'Unit)=0 , (27) 
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where long range interactions are included through the coupling constants K{n — 
m) which are considered to be real, symmetric (i.e. K{n) = K{—n)), and decaying 

oo 

rapidly enough (i.e. such that |i^(m-)| < cxd). As it is clear, in the case of 



m=— oo 



Eq. (HJ it was K{n - m) = 25m,n - Sm,n+l - Sm,n-l- 

The formulated problem (\27\i requires evident modification of the associated 
linear problem and results in a new Green function, which still can be written in 
a form (jB)) but now with the function 

oo 

f{a, t) = e2*<^«-**'^('^), uj{a) = 1 + 2 ^ K{m) cos{2ma) . (28) 

m=l 

Here w(cr) is the dispersion relation of the lattice with long range interactions. The 
new defined function /(a, t) possesses the same properties as the /(a, t) defined 
in © which means that properties ((TB|) and ((T7j) hold also for and thus the 
analogs of the theorems 12. II and |2. 21 can immediately be formulated for the lattice 
with linear long-range interactions. 

Coupling of the next-nearest neighbors can also be nonlinear, which modifies 
the nonlinearity of the Eq. leading to a model 

) = , (29) 

where Wn{-), n € Z, have the properties 

Wn{. . . , Qn-ie"^, Qne"^, Qn+ie"^ , • • • , ^n-ic"**^, g„e"*'^, gn+ie"**^, . . . ) 

= e*'^W„(. . . , qn-l,qn, Qn+l, qn~i,qn, qn+1, ■ ■ ■) , (30) 
Wni^- ■ ■ , Xq-a—m, • • • , Xqn+m, ■ ■ ■ , Xqn—m, ■ ■ ■ , Xqn+m, • • • ) 

— A ^ Wfi(. . . , qji_^, • • • ) qn+m, ■ ■ ■ , qn—m, • • • , qn+m, • • • ) ' (^f ) 

with A and if being an arbitrary real constants, and p being a positive integer 
p>2. 

For the lattice ()29|) . however, not all of the above results allow direct gener- 
alization. Indeed, one essential property used in the preceding sections was the 
conservation of the number of particles, i.e. Qq (or A^, which is the same). But 
not all functions Wn{-) allow this conservation law. More precisely, in order to 
apply the above results of the localization of solution one has to impose one more 
condition which reads 

oo 

[qnWn{. .. ,qn,- ■■ ,qn,- ■■)- qnWn{. .. ,qn, ■■■ ,qn, ■■■)]= . (32) 

n=— oo 

In the same way, to formulate an analog of Theorems 12.11 and 12.21 we have to 
impose one more constraint on Wn{-), namely 

\Wn{...,qn,...,qn,...)\<C\qn\^^P-'^ . (33) 

Then by analogy with the previous considerations we obtain the following 
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Theorem 5.2. Assume that an, n € satisfy the condition moreover that 
\3^) and are valid. Then there exists T > and a unique function q{t) = 
(. . . , . . .), defined in [0,T[ which belongs to the Banach space 

X'^{T), such that qn{t) satisfy 

+00 „t 

qn{t)=Pn{t)+ ^ / G{n - m,t - s)Wn{. ■ ■ ,qn{s), ■ ■ ■ ,q{s)n, ■ ■ ■)ds . (34) 

m=-oo 

This solution can be extended for any T > 0. Such global solution is unique. 

5.3 Saturable nonlinearity 

As the last generalization of the above approach we consider a model with the 
saturable nonlinearity 

ij,.W + A.„W-, ]f2[^^W ^0. (35) 

where A and B are positive constants and x is a real constant. Evidently the 
nonlinear part of this last problem satisfies the conditions and (jSH)) . which 
allows us to formulate 

Theorem 5.3. Assume that an, n £ 1^, satisfy the condition then there exists 
T > and a unique function q{t) = (. . . , g„(t), . . .), defined in 

[0, T[ which belongs to the Banach space X'^{T), such that qn{t) satisfy 

qn{t)=Pn{t) + X £ f G{n-m,t-s)^^^^^^^^ds . (36) 
mt^J^ A + B\qn{s)\^ 

This solution can be extended for any T > 0. Such global solution is unique. 



6 Conclusion 

To conclude, we have proved that initially localized solutions of a number of models 
of the discrete nonlinear Schrodinger type persist localized, with the same degree 
of localization, during an infinite time, and thus constitute a Banach space. The 
proof has been provided using the Banach fixed point theorem, which allowed us 
also to establish the well posedness of the Cauchy problem and thus the impossi- 
bility of existence of collapse. All these properties have rather transparent physical 
explanations which are respectively based on the spectrum of the underlining linear 
lattice having a cut-off frequency, and finite group velocity of the linear excitations 
(in the conventional continuous (nonlinear) Schrodinger equation there is neither 
upper frequency limit nor bound on the speed of the linear excitations). Mean- 
while application of the results was restricted to models which posses an integral 
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of motion in the form of a number of particles N, which guaranteed boundness of 
for all numbers of sites, n, and all times t. Although the most physically 
relevant models do have such an integral of motion, consideration of inter-site non- 
linearities of more general types from the viewpoint of existence and localization 
of the solutions remains an open relevant problem. It is also worth noting that 
other physically important generalizations of the obtained results are related to 
multidimensional and multi-atomic lattices. 
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